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Abstract 

The double differential dilepton spectrum dN/ {dM'j^ dM"^) at fixed transverse 
mass M_|_ allows a direct access to the vector meson spectral functions. Within 
a fireball model the sensitivity of dN/{dMj_ dM"^) against variations of both 
the in-medium properties of mesons and the dynamics of the fireball is investi- 
gated. In contrast to the integrated invariant-mass spectrum dN/dM'^, in the 
spectrum dN/dM'^ dM'j_ with fixed M_|_ the u) signal is clearly seen as bump 
riding on the p background even in case of strong in-medium modifications. 
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I. INTRODUCTION 



A strong motivation for the starting experiments with the HADES detector |I| at the 
heavy-ion synchrotron SIS in GSI Darmstadt is to search for direct evidences of in-medium 
modifications of hadrons via the dilepton (e+e~) channeL Such modifications, hke the p 
meson "melting" and possible in-medium mass shifts, being related to chiral symmetry 
restoration, are considered now as the main cause of the low-mass dilepton excess observed 
in heavy-ion experiments at CERN-SPS energies 0. According to the present understanding 
of the chiral phase transition @,|], the chiral condensate {qq) as order parameter decreases 
with increasing temperature T and baryon density n. This leads to modifications of the 
meson spectral functions since they are coupled to {qq) . While the dependence of the chiral 
condensate as a function of the temperature at small baryon density is quite smooth up to 
a critical temperature of ~ 170 MeV 0, the dropping of {qq) with increasing density at 
small temperature is more pronounced already at moderate densities. For instance, in the 
QCD sum rule analysis 0, which is in accordance with the anticipated Brown- Rho scaling 
0, a density dependence of the hadron mass shift like m* = m(l — 0.18n/no) has been 
predicted, so that even at the nuclear saturation density n = no = 0.15 fm~^ a noticeable 
change of the in-medium mass m* emerges. Unfortunately, in the case of heavy-ion collisions, 
where densities up to Stiq can be expected at SIS energies, the remarkable advantage of the 
strong density dependence meets a problem when trying to identify the predicted mass 
shift by the direct dilepton decay of vector mesons. Due to the collective expansion of the 
baryonic matter the density n depends on time so that obviously a time average of the mass 
shift is to be expected, i.e. the mass shift effect is smeared out over an interval related to 
the change of n. 

To extract the wanted information on the in-medium meson spectral function via dilepton 
measurements in heavy-ion collisions one has also to take into account the general properties 
of the dilepton spectra in a thermalized system which undergoes a collective expansion. For 
instance, in a dense system not too far from local thermal equilibrium, the dilepton rate is 

known to be governed mainly by the Boltzmann factor exp{—E/T}, where E = \J -|- 
is the lepton pair energy, and the pair is characterized by the invariant mass M and three- 
momentum Q. That means the dilepton spectrum dN/dM'^ as a function of the invariant 
mass is a convolution of spectral functions and the Boltzmann factor. For temperatures 
T < 100 MeV, as estimated for SIS energies, and dilepton masses in the region of the 
light vector mesons, M ~ my ~ 770 MeV, the ratio M/T can be in the order of 10 and, 
as a consequence, the resulting spectra are very steep. This can make an identification 
of possible mass shifts, superimposed to (at least collision) broadening, rather difficult. 
Moreover, dynamical effects, caused by flow and cooling, can even stronger mask the wanted 
information. 

In the present note we propose to extend the usual analysis of the integrated spectrum 
dN/dM"^ to the more informative double differential spectrum dN/ {dM]_ dM"^), where M±_ = 
a/ + is the transverse mass of the pair with transverse two-momentum Q±. This 
allows to avoid the majority of the above mentioned difficulties in extracting the in-medium 
vector meson spectral function on the basis of dilepton measurements in heavy-ion collisions. 
In particular, by taking the transverse mass M± at a fixed value (or in a sufficiently narrow 
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bin) one gets rid of the steep slope of the invariant mass spectrum as we show below. As 
a result the spectral function is seen more clearly, and possible in-medium effects become 
more pronounced. 

To demonstrate the very idea of our approach we employ here a fireball model for de- 
scribing the space-time evolution of the dense matter formed in central collisions of heavy 
nuclei at SIS energies. In spite of the schematic character of the model, it is useful to es- 
timate various characteristics, such as the relative contributions of vector mesons decaying 
inside and outside of the fireball and Dalitz decays of mesons and the A as well, which 
are related to observables accessible in HADES experiments. Along with more complicated 
transport model calculations we consider our approach as a first step towards elucidat- 
ing the most appropriate observables for verifying in-medium changes of meson properties 
in heavy-ion collisions. 

Our paper is organized as follows. In Sect. 2, we present the parameterizations of the 
dilepton rates which we use in our exploratory studies. The fireball dynamics is given in 
Sect. 3. The resulting dilepton spectra are discussed in Sect. 4, where we contrast the 
dileptons emitted during the fireball's lifetime and such dileptons which stem from hadron 
decays after freeze-out. The discussion can be found in sect. 5 together with remarks on the 
support of the presented results by calculations based on a BUU code. 



II. DILEPTON RATES 



The dilepton rate from vector meson decays, V —>■ e^e with V = p,u},(j), in an ideal 
resonance gas approximation is given by 

dN ^ 2dv ( u-Q\l (myr*?*)(Mr 



^ (27r)3 """"P \ T j TT (M2 - ml^y + (myP/)^ ' ^ ^ 



d^xd^Q Y^27r)3 T j rr {M^ - m^y)^ + {mvT 

where Fy* and Ty^^+e- denote the total and dilepton decay widths, my is the vector meson 
mass, and dy stands for the degeneracy factor. The lepton pair has the four-momentum 
Qf^ = (M_|_ coshy, sinhy, Q^), where Y denotes the pair rapidity. In (|I|) we use the 
Breit-Wigner parameterization of the meson spectral function. 



1 myT 



tot 
V 



^^(^) = r .^2_^2,2^r^,.rtot,2 > (2) 



TT (m^ — my)2 + {myT 



V 



which can be thought to result from the generalization of an on-mass shell delta function to 
the case of an unstable particle [pH]. Such an ansatz is also used for QCD sum rule analyzes 



111 of vector meson properties in a nuclear medium. Due to four-momentum conservation 
one gets immediately = Q"^ = M^. 

It should be emphasized that, within the vector dominance model, the dependences 
rp_+e+e- oc and r*°* oc just cancel, therefore, in (|I|) the widths are understood to be 

constant. A remaining M dependence stemming from r*'^^+^_ oc (1 — )^^^ is partially 
mimicked by the factor M in the numerator in . We have checked that the rate described 
by (|l|) agrees with the standard parameterization ||10| of the process vr+vr" —>■ e^e~ in the 
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intervals 400 MeVQ < M < 1000 MeV and 50 MeV < T < 100 MeV. For the narrow u and 
mesons these subtleties are not important. 

The collective expansion of the matter is characterized by the four-velocity u^{x) de- 
pending generally on space and time, as also does the temperature T(x). In spite of the 
approximative character of (|T]) it is useful for demonstrating how the in-medium modifica- 
tions in the spectral function @ show up in the dilepton spectra if flow effects of the matter 
are taken into account. 

Integrating the rate (|l]) over the three-momentum of the lepton pair one gets the rate as 
a function of the invariant mass, 

dN ^ dy (M\ (mv.rV°^)(Mr^^,+e~) . . 

d^x dM^ ^ 27r3 ) {M^ - + (my 1*^)2 ' ^ > 

where Ki is a Bessel functions. 

In Fig. la we plot the rate @ at a temperature of T = 70 MeV to demonstrate the 
general tendency of the very steeply dropping spectrum with increasing M. Since we do not 
advocate in the present note any particular functional form of the meson spectral function 
(2), which is under intense debate so far, we parameterize schematically here the in- medium 
width broadening by the change of the total vacuum width according to Fy — > Fy = 6Fy 
with a constant parameter b = 3 for all vector mesons. (The position of the in-medium 
vector meson peak is not important for this example.) As seen in Fig. la the important 
effect of the vector meson broadening is smeared out by the rapidly dropping thermal weight 
factor Ki{M/T) in (0). No peak related to the p meson remains. The u peak is still visible, 
but not very pronounced. The strongest peak structure is related to the meson. 

The above example can be contrasted to the double differential rate dN/ {d^xdM]_dM'^) 
obtained from (0) by integrating over the dilepton rapidity and the polar angle. For a 
spherically symmetric source we find 

dN ^ 2dv (M^^\ f Q^Vrl \ (myFV°^)(MFy^,+e-) ... 

d'^xdMldM^ V(2vr)3 °V ^ ° V T ) {M^ - mlf + {myV^^Y' 

where 7 = 1/ a/1 — is the Lorentz factor corresponding to the radial flow three- velocity 
Vr, and = M\ — . One can infer from (^ that, by fixing the transverse mass, the 
thermal weight is frozen in. For sufficiently small values of the flow velocity the rate 
dN I {d'^xdM]_ dM"^) as function of the invariant mass at fixed M_|_ reflects directly the shape 
of the meson spectral function. In Fig. lb we show this remarkable property for Mj_ = 1.3 
GeV, Vr = 0.3 and T = 70 MeV. The p bump is here clearly visible below the u peak; the 
(j) needle is more pronounced. The price is, of course, a diminished count rate. However, for 



At this and smaller values of M the rate is anyhow dominated by Dalitz decays. Therefore, we 
do not attempt a modification of the low-M tail of (|l|) which arises from a particular model. For 
more involved expressions of the rate we refer the interested reader, e.g., to ||l2|, where also the 
coupling to the baryon charge is included. 
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M±_ = 0.9 GeV the rate is still not too much smaller in comparison with the rate displayed 
in in Fig. la at M = 0.9 GeV. 

In Fig. lb we show that the rate increases in the p region with increasing values of the 
radial velocity Vr- This fact is known since some time [0 as enhancement of the p peak 
due to transverse flow. Note that the rate dN / {d'^xdM'^) is completely unaffected by flow 
effects. 

It should be stressed that other representations of the double differential rate, such as 
dN/d'^x dQ\dM'^ , are less suitable for an access to the spectral function. 



III. DYNAMICS 

To obtain the dilepton spectrum from the meson decay rate (|l|) one needs to specify the 
space-time evolution of the matter. According to our experience Jl^ the detailed knowledge 
of the space-time evolution is not necessary to describe the experimental dilepton spectra in 
relativistic heavy-ion collisions in the low-mass region, where the in-medium modifications 
of the meson spectral function play a crucial role. In this line, and also to avoid too many 
parameters, we employ here a variant of the blast wave model |T3 with constant radial 



expansion velocity Vr- In such a model the radius of the fireball increases like R{t) = 
Ri + Vrt. From the analysis of hadron ratios in heavy-ion reactions at SIS energies a freeze- 
out temperature of Tf = 50 MeV at freeze-out baryon density n/ = 0.3no is obtained 
||16|| . Extrapolating back along an isentropic line by using a resonance gas model as in 
||T6| up to rii = 3no (as suggested by transport code simulations for the maximum density 
at SIS energies) one finds Tj = 90 MeV. The freeze-out time t/ is determined by baryon 
conservation within the fireball volume; we use as baryon participant number = 330. 
The time evolution starts at initial time ti = 0. We take Vr = 0.3 as a reasonable average 
according to |M. By these parameters the dynamics of the fireball is completely fixed. 



IV. DILEPTON SPECTRA 
A. The fireball radiation 

Since the aim of our work is a demonstration of the advantage of the double differential 
dilepton spectrum dN/{dM\dM'^) independently of a given form of the vector meson's 
spectral function Ay, we do not use specific functional changes of Ay but consider a variety 
of possible parameterizations of the in-medium effects via mass shifts and width broadenings, 
i.e., changes of the parameters Fy* and my. As long as the quantitative description of Ay{m) 
in a medium is matter of debate P] such a procedure seems appropriate. 

We begin with an extreme case where the p meson mass shift is Amp = 300 MeV, 
while the u meson is unaffected, i.e. Am^ = 0, and my = my — Amy. The in-medium 
widths of both the low-lying vector mesons is increased by 6 = 3. The resulting spectra are 
displayed in Fig. 2, and we have chosen M±_ = 0.9 GeV. In contrast to the integrated rate 
dN/dM"^, where the p peak disappears (cf. Fig. 2a), one can clearly observe the splitting of 
both vector meson peaks in the dN/{dM]_dM'^) spectrum (cf. Fig. 2b), thus reflecting the 
original modifications of the p, uj spectral functions. 
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Another similarly extreme case is to keep the p mass fixed and to shift the uo mass by 
Am^^ = 300 MeV down; again with widths broadening factor 6 = 3. As seen in Fig. 3 the 
information on the change of the meson properties is best visible in the dN/ {dM]_ dM'^) 
spectrum, where a strong u peak and some part of the p bump are distinguishable (cf. 
Fig 3b). Otherwise, the shifted u peak is also nicely visible in the dN/dM"^ spectrum (cf. 
Fig. 3a) because it "rides" on the p background and becomes more pronounced; a p bump 
is here again hardly visible. 

In more realistic considerations the values of the meson mass shifts and width broad- 
enings should depend on the baryon density and temperature which are strongly time de- 
pendent. This in turn causes an additional smearing of the peaks of the spectral functions. 
To demonstrate that dN/{dM]_dM'^) is still useful to elucidate occurrence of in-medium 
modifications we parameterize now the in-medium u mass by m* = muj{l — 6^ ■ n/no), 
-ptot* _ -ptot _|_ ^Y^^n/riQ with 6^ = 0.08 and ^F^^ = 20 MeV, and for the p meson we assume 
only a broadening according to F*°** = F*°* + STp ■ n/nQ with STp = 400 MeV, where the 
superscript * denotes again the in-medium quantities. These parameterizations are in line 
with current expectations ||T^ that the p meson suffers a very strong in-medium broadening, 
while for the u meson one can expect moderate modifications. The resulting spectra are dis- 
played in Fig. 4. In contrast to the integrated invariant mass spectrum dN/dM"^ (cf. Fig. 4a), 
in the spectrum dN/{dM]_dM'^) with fixed M± (cf. Fig. 4b) the u meson contribution is 
clearly seen as shifted bump riding on the p background, in spite of the smearing caused by 
the space time evolution. (Notice that in Fig. 4a the smeared u contribution could easily be 
misinterpreted as p bump.) Therefore, the measurement of the double differential dilepton 
spectrum at fixed M± offers much better chances to identify a possible u mass shift. 

It should be stressed that changes of Tj cause a noticeable up/down shift of the spectra 
but let the shapes nearly the same. 



B. Background contributions 



Within the given model one can also estimate various background contributions. For 
instance, the contributions of vector meson decays after freeze-out can be consistently eval- 
uated. The invariant mass spectrum of such decays can be calculated from 



dN 
dAP 



1 dr 



V- 



/ J ptot 



dM2 



-Ny{T^,my = M), 



(5) 



-/dM"^ = Tv^e+e-^viM'^) and Ny{Tf,my) is the number of vector mesons 
This number is calculable within the Cooper- Frye formalism |]18[, which 



where dVy^e+i 
at freeze-out. 

includes relativistic effects related to the matter's expansion. In our spherical fireball model 
with constant radial fiow velocity one gets 



Ny(Tf,my 



2 ^2 + (myrft^2 




(6) 



+ Tf) — Tf cosh a/ 
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where m± = a/ m'^ + is the transverse mass of the meson V with transverse momentum 
p±, and df = 'yVrP±/Tf. In deriving (^) we employ the meson spectral function Av{m) from 

(I)- 

The calculations show that even for the "long-living" u, (j) mesons the decay contribution 
after freeze-out is considerably smaller than the decay contribution during the fireball's 
life time. This is in contrast to naive expectations based on simple life time arguments, 
according to which the life time of the uj meson is 23 fm/c compared to the life time of 
the thermalized system (9(10) fm/c. To arrive at a consistent estimate one has to take into 
account the evolution and the hypothesis of thermalization. The latter assumption is the 
key to understand the small contribution of the decays after freeze-out: In the early stage 
the temperature is much higher and, according to the assumed chemical equilibrium, the 
appearance of vector mesons is much more frequent than at freeze-out. At the pole positions 
the ratio of the fireball to the freeze-out decay contributions can be roughly estimated by 
exp{— my/Tj}/ exp{— my/Tj}, because of the dominating exponential factors. 

It should be emphasized that the analysis of the hadron multiplicities at freeze-out in 
|T6| supports the chemical equilibrium picture for some hadrons. In the present work we 
assume chemical equilibrium for the full time evolution. This equilibrium is certainly not 
maintained during the time evolution for all hadron species. For instance, the uj production 
via the inverse decay reaction Svr — is rather unlikely. Such effects can be simulated by 
additional chemical off-equilibrium potentials, which keep certain hadron species on a higher 
number than we have assumed here. In this respect our estimates, e.g., of the uj numbers, 
serve as unfavorable lower limits when assuming particle production in the stages of highest 
temperature and density. 

While we focus on the direct decay channels V —>■ e~^e~ , the Dalitz decay channels, like 
A Ne~^e~, (p — > r]e~^e~, uj — > Tr^e~^e~, rj' — > Tr'^{r])e~^e~ , r] — > •ye'^e', ttq 7e"'"e~ etc., 
populate the invariant mass spectrum in the low-mass region (for a recent comprehensive 
analysis of the various other channels cf. |]19|)- These Dalitz decays are calculated by (5, 6) 
with dTy^e+e- /dM"^ to be replaced by dTA-.Be+e-fdM'^ (for hadron decays after freeze-out) 



dNA^Be+e- ^A f , 2 mAF*/ dT A^Be+e- 



and 



d^xdM^ 27r3 J {m^ - m\f + {rriAV^l^Y dM^ 

dTA^Be+e- _ Q -p f^s 

dM^ "3^^^^^^*' 

(for hadron decays during the fireball's life time ), where Ta^by* is taken from for meson 
decays and from for A decays. 

With respect to the A Dalitz contribution we mention that the number of A resonances 
depends on the baryo-chemical potential. The latter one is in our fireball model approx- 
imately 850 MeV and fairly well time independent. The quoted value results from the 
chemical freeze-out analysis in [0, and the approximate time independence is a special 
feature of the isentropic line with specific entropy of 5 in the interval n = 0.3 ■ ■ ■ Suq. 

We have checked that the sum of these Dalitz decay channels at M > 0.5 GeV is below 
the direct vector meson decay contributions (see also remarks in the next section). 
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V. DISCUSSION AND SUMMARY 



In summary we have demonstrated that, within a fireball model and a simple parame- 
terization of the rate, the double differential dilepton spectra dN / {dMj_ dM"^) at fixed value 
of M_[_ can deliver valuable information on a possible change of the in-medium properties of 
vector mesons. The usually considered invariant mass spectrum dN/dM"^ contains a convo- 
lution of the spectral function with a thermal weighting function which is too strong such 
that modifications of the spectral functions are difficult to disentangle, in particular if the 
dynamics of the system causes an additional smearing. 

We consider our fireball model as first step towards finding the most appropriate char- 
acteristics of dilepton spectra to probe in-medium vector meson spectral functions. The 
results obtained, in particular the dN/dM"^ dM]_ spectra, should be confirmed by more real- 
istic transport calculations. Indeed, using the BUU code developed from ||20| we performed 
calculations in line with the present work. These calculations show very clearly that, when 
selecting the interval 850 MeV < < 950 MeV, (i) the resulting spectrum dN/dM'^ dMi 
is dominated by the channel tt^tt^ ^ p ^ e^e~ in the interval 300 MeV < M < 850 MeV 
thus directly reflecting the shape of the p formfactor, (ii) in the high-M tail also direct p 
decays, with p's stemming from other channels than vr+vr" annihilations, can become as 
large as the 7r"'"7r~ annihilation contribution, (iii) the uj rides on the p bump, and (iv) the pn 
bremsstrahlung, A Dalitz and rj Dalitz decays are negligible in the mentioned M interval. 
The details of these BUU calculations will be presented elsewhere . 

To make reliable predictions for the future experiments at HADES one needs also addi- 
tional analyses of the mass shifts and widths with respect to both the baryon density and 
temperature dependences in the regions expected to be achieved in heavy-ion collisions at 
SIS. 
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FIGURES 




M[GeV] M[GeV] 



FIG. 1. The rates dN / {d^x dW^) ((a) left panel) and dN / {d^x dMldM"^) ((b) right panel, 
M_L = 1.3 GeV) for T = 70 MeV, Vr = 0.1 and 0.3, and 6 = 3. 




M[GeV] M[GeV] 

FIG. 2. The rates dN/{dM^) ((a) left panel) and dN/{dMl dM"^) ((b) right panel, M_l = 0.9 
GeV) for the dynamical scenario described in the text and for 6 = 3, Amp = 300 MeV, Am^^ = 0. 
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FIG. 3. As in Fig. | but for Arrip = 0, Am^ = 300 MeV. 




M[GeV] M[GeV] 

FIG. 4. As in Fig. | but for 6^ = 0.08, 5T^ = 20 MeV, 6p = 0, 6Tp = 400 MeV. The uj (p) 
contribution is displayed separately as short-dashed (long-dashed) line. 
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